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ABSTRACT 

For  life  assessment  of  airframe  components,  it  is  important  to  understand  how  the  residual 
stresses  due  to  material  yielding  affect  crack  initiation  or  crack  growth  at  holes.  A  key  step  in 
understanding  such  effects  is  the  quantification  of  the  Mode  1  stress  intensity  factor  for  cracks 
that  may  initiate  subsequent  to  overloading.  Here,  a  two-dimensional  weight  function 
approach  is  used  to  determine  stress  intensity  factors  for  cracks  in  either  tensile  or 
compressive  stress  fields  due  to  one  of  three  mechanisms:  remote  tension  overload,  remote 
compression  overload  or  hole  cold  expansion.  The  effect  of  subsequent  remote  loading  is  also 
considered.  The  key  input  is  the  stress  distribution  in  the  corresponding  uncracked  body 
along  the  prospective  crack  path.  The  key  trends  are  investigated  through  many  numerical 
examples  for  symmetrically-cracked  holes  in  large  steel  and /  or  aluminium  alloy  plates.  Cold 
expansion  of  finite-width  plates  representative  of  C-130  wing  skin  locations  is  also  studied. 
For  both  remote  overload  cases,  it  is  shown  that,  once  the  crack  length  is  the  same  or  larger 
than  the  initial  yield  zone,  the  stress  intensity  factors  are  the  same  as  for  the  case  without  the 
initial  overload.  However,  for  cold-expanded  holes,  the  beneficial  reduction  in  stress  intensity 
factor  extends  beyond  the  initial  yield  zone.  Hence,  the  present  work  provides  a  greater  depth 
of  understanding  of  how  typical  residual  stresses  can  affect  key  inputs  into  airframe  life 
assessment. 
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Investigation  of  Stress  Intensity  Factor  for 
Overloaded  Holes  and  Cold-Expanded  Holes 

Executive  Summary 


Significant  work  undertaken  in  Air  Vehicles  Division  involves  providing  advice  on 
through-life  support  of  Australian  Defence  Force  aircraft.  For  accurate  life  assessment  of 
airframe  components,  it  is  important  to  have  an  understanding  of  how  the  residual 
stresses  due  to  material  yielding  affect  crack  initiation  or  crack  growth  at  holes  in  metallic 
components.  A  key  step  in  understanding  such  effects  is  the  quantification  of  the  Mode  1 
stress  intensity  factor,  for  cracks  which  may  initiate  subsequent  to  the  overloading.  The 
weight  function  method  is  a  useful  approach  to  determine  Mode  1  stress  intensity  factors, 
since  the  crack  does  not  need  to  be  explicitly  modelled.  Hence,  in  the  present  investigation, 
a  two-dimensional  weight  function  approach  is  used  to  determine  stress  intensity  factors 
for  cracks  in  either  tensile  or  compressive  stress  fields,  due  to  one  of  three  mechanisms: 
remote  tension  overload,  remote  compression  overload  or  hole  cold  expansion.  An 
important  analysis  input  is  the  stress  distribution  in  the  corresponding  uncracked  body 
along  the  prospective  crack  path,  obtained  from  either  finite  element  analysis  or 
theoretical  methods. 

In  the  present  work,  analyses  are  undertaken  to  consider  key  trends  and  limitations  of 
current  prediction  methods,  with  an  emphasis  on  the  use  of  weight  function  methods.  The 
main  areas  addressed  are  as  follows.  Initially,  a  detailed  literature  survey  is  presented  on 
the  use  of  weight  function  approaches  in  the  context  of  three  key  types  of  overloads: 
remote  tensile  overload,  remote  compressive  overload  and  local  cold  expansion.  As  an 
example  for  the  generation  of  stresses  along  the  prospective  crack  path,  finite  element 
analysis  of  an  uncracked  hole  in  a  remotely-loaded  steel  plate  is  given.  Stress  intensity 
factor  results  for  three  sets  of  numerical  test  cases  are  then  determined.  Two  of  these 
involve  cold-expanded  holes  in  large  D6ac  steel  plates  and  large  aluminium  alloy  plates. 
The  third  test  case  geometry  corresponds  to  a  fatigue  test  coupon  consisting  of  a  finite- 
width  aluminium  alloy  plate,  which  is  representative  of  locations  on  the  lower  wing  skin 
of  the  C-130  aircraft.  Finally,  results  and  comparison  of  generic  cases  for  a  D6ac  steel  plate 
for  all  three  overload  cases  are  given. 

For  both  remote  overload  cases,  it  is  shown  that,  once  the  crack  length  is  the  same  or  larger 
than  the  initial  yield  zone,  the  stress  intensity  factors  are  the  same  as  for  the  case  without 
the  initial  overload.  However,  for  cold-expanded  holes,  the  beneficial  reduction  in  stress 
intensity  factor  extends  beyond  the  initial  yield  zone.  A  set  of  Fortran  computer  programs 
has  been  written  to  generate  the  results  in  this  report.  These  programs  are  given  on  an 
enclosed  Compact  Disc  for  future  reference  and  use.  The  present  work  has  provided  a 
greater  depth  of  understanding  of  how  typical  residual  stresses  can  affect  key  inputs  into 
airframe  life  assessment.  The  use  of  these  techniques  can  assist  in  providing  more  timely 
and  accurate  prediction  of  fatigue  life,  for  Australian  Defence  Force  platforms. 
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1.  Introduction 

Holes  are  the  most  common  stress  concentration  feature  in  metallic  aircraft  structures;  hence 
the  growth  of  fatigue  cracks  at  holes  is  a  common  occurrence.  Overloads  at  holes  can  occur 
due  to  large  service  loads  experienced  by  airframe  components.  These  loads  can  be  tensile  or 
compressive,  and  may  result  in  material  yielding  at  the  hole  edge.  Compressive  overloads 
may  more  commonly  occur  on  upper  wing  surfaces  while  tensile  overloads  typically  occur  on 
lower  wing  surfaces.  A  tension  overload  causes  local  yielding,  followed  by  unloading,  can 
provide  beneficial  compressive  residual  stresses  at  the  hole  edge.  Conversely,  compressive 
overloads  followed  by  unloading  can  lead  to  potentially  detrimental  tensile  residual  stresses. 
Holes  can  also  be  overloaded  locally  to  enhance  fatigue  life.  For  example,  cold  working  of 
holes  is  a  technique  that  is  widely  used  in  aircraft  structures  to  induce  beneficial  compressive 
residual  stresses  around  the  hole.  This  is  achieved  using  a  mandrel  drawn  through  the 
starting  hole  in  order  to  apply  radial  pressure  to  the  hole  surface,  causing  local  yielding  and 
material  elongation  in  the  circumferential  direction. 

For  the  life  assessment  of  airframe  components,  an  understanding  of  how  the  residual  stresses 
and  initial  yield  zone  affect  crack  initiation  or  crack  growth  is  important.  A  typical  step  in 
understanding  such  effects  is  the  quantification  of  the  Mode  1  stress  intensity  factor,  if, ,  for 
cracks  which  may  initiate  subsequent  to  the  overloading.  It  is  desirable  to  have  available  some 
knowledge,  even  trends,  of  the  following:  (i)  conditions  for  which  overload  effects  are 
significant;  (ii)  how  significant  the  effects  are  for  different  cases,  for  example,  how  the  effects 
due  to  remote  overloads  compare  to  local  treatments  such  as  cold  expansion;  and  (iii)  the 
accuracy  and /  or  limitations  of  common  analysis  approaches  for  predicting  K t .  Based  on  a 
review  of  the  literature,  it  appears  that  only  limited  work  has  been  conducted  to  address  these 
issues.  It  is  also  clear  that  most  relevant  work  has  focussed  on  cold  expansion  of  holes,  and  on 
the  application  of  weight  function  methods  to  holes  or  other  crack  configurations. 

Hence,  in  the  present  work,  analyses  are  undertaken  to  consider  key  trends  and  limitations  of 
current  prediction  methods,  with  an  emphasis  on  the  use  of  weight  function  methods.  A 
detailed  literature  survey  is  presented  in  Section  2,  covering  the  use  of  weight  function 
approaches  in  the  context  of  three  key  types  of  overloads:  remote  tensile  overload,  remote 
compressive  overload  and  local  cold  expansion.  The  general  plate  geometry  and  loading 
arrangement  used  in  the  computational  analyses  are  given  in  Section  3.  The  finite  element 
analysis  (FEA)  of  an  uncracked  hole  in  a  remotely  loaded  plate  is  presented  in  Section  4.  In 
Section  5,  the  weight  function  approach  is  given  for  cases  of  a  hole  in  a  remotely-loaded  plate, 
as  well  as  for  a  cold-expanded  hole.  Section  6  then  gives  the  stress  intensity  factor  results  for 
three  sets  of  numerical  test  cases.  Two  of  these  involve  cold-expanded  holes  in  large  D6ac 
steel  plates  and  large  aluminium  alloy  plates.  The  third  test  case  geometry  corresponds  to  a 
fatigue  test  coupon  consisting  of  a  finite- width  aluminium  alloy  plate,  which  is  representative 
of  locations  on  the  lower  wing  skin  of  the  C-130  aircraft.  Then  Sections  7,  8  and  9  give  the 
results  and  comparison  of  generic  cases  for  a  D6ac  steel  plate  for  all  three  overload  cases. 
Here,  Section  9  focusses  on  a  comparison  of  results  for  the  cold-expanded  hole  and  the 
tension-overload  case.  The  conclusions  are  then  given  in  Section  10.  A  Compact  Disc 
containing  source  code  for  the  numerical  integration  and  the  examples  is  provided  with  this 
report. 
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2.  Literature  survey 

Weight  function  techniques  for  estimating  K f  were  first  conceived  by  Bueckner  [1],  Rice  [2] 
also  obtained  the  same  equation  by  a  different  method  and  demonstrated  that  the  weight 
function  is  a  universal  function.  The  method  is  very  general,  since  estimating  Kl  depends  only 
on  the  integration  of  the  weight  function  and  the  stress  distribution  along  the  prospective 
crack  path  in  the  corresponding  uncracked  body. 


2.1  Weight  functions  and  their  application 

A  large  number  of  weight  functions  have  been  determined  for  different  conditions,  such  as 
cracks  emanating  from  holes,  cracks  subject  to  non-uniform  stress  fields  and  cold  working  of 
holes  containing  cracks  (see  for  example  Wu  and  Carlsson  [3]).  Some  further  examples  are  as 
follows.  Kiciak  et  al  [4]  have  derived  weight  functions  for  corner-elliptical  cracks  in  a  finite¬ 
thickness  plate.  Chen  and  Wang  [5]  have  derived  weight  functions  for  cracks  subject  to  non- 
uniform  stress  fields  and  cold  working  of  holes  containing  cracks.  Glinka  and  Shen  [6]  have 
derived  weight  functions  for  edge  cracks  in  infinite  plates  and  through  cracks  in  infinite 
plates.  Mawatari  and  Nelson  [7]  have  used  Gauss-Chebyshev  integration  in  order  to  eliminate 
the  singularity  and  reduce  the  error  below  1  % .  Similarly,  Daniewicz  [8]  has  used  the  same 
method  and  has  achieved  increased  computational  efficiency.  Grandt  and  Kullgren  [9]  have 
tabulated  stress  intensity  factors  for  flawed  fastener  holes  using  a  weight  function  technique. 
Paris  et  al  [10]  have  derived  weight  functions  for  a  number  of  two-dimensional  (2D)  problems. 
Shen  and  Glinka  [11]  have  derived  weight  functions  for  2D  surface  semi-elliptical  cracks  in 
three-dimensional  (3D)  plates.  Petroski  and  Achenbach  [12]  have  computed  the  weight 
function  from  stress  intensity  factors  for  a  reference  loading  configuration.  They  obtained 
crack  face  displacements  that  were  in  good  agreement  with  analytical  results  for  a  number  of 
cases,  and  their  technique  involves  simple  quadrature  and  is  very  efficient.  Sobczyk  and 
Trqbicki  [13]  looked  at  the  effect  of  random  residual  stresses  on  fatigue  crack  growth  and  used 
the  weight  function  technique.  A  stochastic  fatigue  crack  growth  equation  was  also  used  to 
determine  the  random  number  of  cycles  to  failure.  Jones  et  al  [14]  have  proposed  weight 
functions  for  2D  and  semi-elliptical  cracks  emanating  from  a  notch,  which  are  applicable  for 
both  short  and  long  cracks.  Zhao  and  Wu  [15]  have  evaluated  stress  intensity  factors  for  a 
semi-elliptical  surface  crack  in  a  semi-circular  notch  by  use  of  the  weight  function  technique 
for  loads  applied  on  the  crack  face  and  the  edge  of  the  specimen.  Good  results  are  obtained  in 
comparison  to  FEA.  Work  undertaken  by  Roy  and  Saha  [16]  involved  the  development  of  an 
analytical  expression  for  the  weight  function  of  an  elliptical  crack  under  normal  loading,  and 
the  solution  is  verified  for  three  separate  cases. 


2.2  Cold  expansion 

For  the  case  of  plates  with  cold-worked  holes  subject  to  remote  loading,  Pinho  et  al  [17]  have 
used  2D  analytical  and  semi-analytical  methods  to  predict  stress  intensity  factors.  Here  a  focus 
was  the  study  of  loading  conditions  and  yield  zone  size  required  to  open  the  crack  fully  for  Al 
2024  Alclad  aluminium  alloy  material.  Clark  [18]  has  used  the  weight  function  technique  to 
calculate  residual  stress.  Here  comparisons  with  experimental  results  have  shown  that  the 
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method  can  give  reasonable  predictions  of  critical  crack  length.  Pavier  et  al  [19]  have  used  2D 
and  3D  FEA  and  2D  weight  functions  to  investigate  cold  expansion  of  cracked  holes. 
Normalised  results  are  given  showing  that  the  weight  function  technique  underestimates  the 
stress  intensity  factor.  The  agreement  improves  with  increasing  crack  length  and  when  the 
remote  loads  are  high  compared  to  the  material  yield  point.  Moreira  et  al  [20]  have  used  both 
FEA  and  the  weight  function  technique  to  determine  residual  stress  intensity  for  cracked  cold- 
expanded  holes  (i.e.  with  no  subsequent  remote  loading)  under  plane  stress  conditions.  Here 
good  agreement  between  both  methods  for  a  range  of  crack  sizes  was  shown. 


2.3  Compressive  overloads 

Very  few  references  exist  on  the  effects  of  compressive  overloads  applied  to  plates  with  open 
holes.  However,  heat  treatment  on  1080  steel  has  been  used  to  induce  tensile  residual  stresses 
in  a  notch  in  a  compact  tension  specimen  by  Prawoto  [21].  The  predictions  of  stress  intensity 
based  on  linear  elastic  fracture  mechanics  and  weight  functions  compared  well  with  the 
experimental  results. 


2.4  Tensile  overloads 

In  the  case  of  tension  overloads  for  plates  with  open  holes,  some  references  do  exist.  An 
experimental  study  has  been  carried  out  by  Wahab  et  al  [22],  It  was  found  that  tensile 
overloads  retard  the  crack  growth  rate  for  a  period  of  time.  However,  when  the  crack  grows 
out  of  the  plastic  zone  the  crack  growth  is  no  longer  retarded.  Similarly,  the  residual  stresses 
considered  by  Wang  [23]  are  tensile  in  the  plate,  and  fatigue  crack  growth  is  promoted  when 
the  crack  grows  into  the  residual  stress  field.  Also,  for  the  tension  overload  case  LaRue  and 
Daniewicz  [24]  have  used  an  elastic-plastic  FEA  to  relatively  accurately  predict  fatigue  crack 
growth  using  a  weight  function  computation.  From  references  [19, 20, 21],  there  appears  to  be 
some  inconsistency  in  the  claimed  effectiveness  of  weight  function  methods  to  predict  the 
stress  intensity  factor  when  the  crack  is  in  a  local  compressive  residual  stress  field.  Work 
carried  out  by  Yu  and  Abel  [25]  has  developed  a  weight  function  technique  capable  of 
predicting  stress  intensity  during  partial  crack  closure.  This  involves  computing  the  crack 
surface  displacement  and  a  trial-and-error  method  to  calculate  K f ,  and  is  applied  to  a  centre- 
cracked  plate.  The  method  is  fairly  general  and  may  be  an  option  for  application  to  cracking  at 
holes. 


3.  Geometry  and  loading  arrangement 

In  this  section,  the  geometry  and  notation  for  the  two  general  cases  considered  are  given  as 
follows:  (i)  overload  due  to  remote  loading  of  a  plate  containing  a  central  hole  and  (ii) 
overload  due  to  cold  expansion  of  a  centrally-located  hole  in  a  plate.  In  both  cases  subsequent 
elastic  remote  loading  can  be  introduced. 
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3.1  Remotely-loaded  plate  with  hole 

A 1  mm  thick  plate  of  width  2 w  and  height  2 h  containing  a  central  hole  of  radius  R  is  shown 
in  Figure  1.  The  plate  is  subject  to  an  initial  remote  load  <jr  which  may  be  compressive  or 
tensile.  Unless  noted  otherwise,  all  analyses  are  for  a  large  plate,  typically  with  w  =  1  OR. 


<?r 


Figure  1:  Plate  with  central  hole  loaded  by  remote  stress 


3.2  Cold-expanded  hole 

The  cold-expanded  hole  is  idealised  as  shown  Figure  2.  Once  again,  the  plate  width  is  given  as 
w,  and  here  the  plate  is  considered  infinitely  long.  The  interference  level  is  defined  by: 


X  -  Rm  ~  R 

R 


(1) 


where  RM  is  the  initial  radius  of  the  mandrel  and  R  is  the  initial  radius  of  the  hole. 
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r 


°r 

Figure  2:  Cold-expanded  hole  in  an  infinitely -long  plate  of  finite  width  subjected  to  remote  stress 

4.  FEA  for  uncracked  hole  in  remotely-loaded  steel  plate 

The  following  sub-sections  cover  the  finite  element  mesh,  modelling  parameters,  boundary 
conditions  and  typical  stress  results  due  to  compressive  and  tensile  overloads.  Some  results 
corresponding  to  the  subsequent  application  of  a  tensile  service  load  are  also  given.  The 
analysis  is  two-dimensional  with  plane  stress  conditions  being  assumed  unless  otherwise 
stated.  The  ABAQUS  version  6.9  structural  analysis  software  code  is  used  for  the  FEA.  Unless 
noted  otherwise,  all  analyses  represent  a  large  D6ac  steel  plate.  These  stress  results  are  then 
used  in  Sections  7  and  8  to  determine  stress  intensity  factors  using  the  weight  function 
method  for  tensile  and  compressive  overload  cases  respectively. 
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4.1  Finite  element  mesh  and  material  modelling  parameters 

4.1.1  Mesh  and  elements 

Consider  a  uniaxially-loaded  square  plate  with  a  central  hole  as  shown  in  Figure  3,  which 
leads  to  the  development  of  regions  of  plastic  flow  as  shown  in  Figure  4.  The  finite  element 
model  assumes  Vi  symmetry.  As  shown  in  Figures  5  and  6,  the  region  along  the  prospective 
crack  path  has  a  highly-refined  mesh  consisting  of  4-noded  quadrilateral  elements.  Each 
element  has  an  edge  length  of  0.01  mm  compared  with  the  hole  diameter  of  10  mm.  Three- 
noded  triangular  elements  have  been  used  for  the  transitioning  of  the  mesh  and  the  remainder 
of  the  Vi  model. 


ar 

Figure  3:  Geometry  for  finite  element  model  showing  prospective  crack  location 
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Crack  length  controlled  with  symmetry  restraints 


4.1.2  Material 

The  analysis  is  materially  nonlinear  and  the  von  Mises  yield  criterion  is  used.  The  material 
used  is  D6ac  steel  and  the  properties  are  given  in  Table  1  (see  [26]). 

Table  1:  Material  properties  for  constitutive  model  for  D6ac  steel 


E  (MPa) 

o>  (MPa) 

EP  (MPa) 

200,000 

1,309 

1,824 

In  the  above  Table,  E  is  the  Young's  Modulus,  <jy  is  the  yield  stress,  and  Ep  is  the  modulus 
in  the  plastic  range.  It  can  be  seen  that  the  response  is  close  to  elastic-perfectly  plastic. 

4.1.3  Boundary  conditions 

Displacement  constraints  due  to  symmetry  are  applied  in  the  A  direction  along  the  left-hand 
boundary  at  X=  —R,  as  shown  in  Figure  5.  Similarly,  displacement  constraints  are  also  applied 
in  the  Y  direction  along  the  top  boundary  at  Y=  0  in  Figure  5. 

4.1.4  Loading  and  nonlinear  convergence 

The  ABAQUS  code  automatically  selects  the  number  of  load  steps  to  be  used  in  the  analysis  to 
achieve  a  selected  convergence  criterion.  Here,  for  an  increment  of  load,  the  sum  of  the 
internal  nodal  forces  is  compared  to  the  total  applied  load.  The  criterion  chosen  is  that  for  an 
acceptable  solution  this  difference  is  less  than  0.5%. 

4.1.5  Mesh  convergence  check 

Since  the  FEA  mesh  is  very  fine,  no  convergence  testing  based  on  mesh  refinement  was 
carried  out.  However,  certain  accuracy  checks  have  been  undertaken.  For  example,  for  the 
elastic  model  the  maximum  discontinuity  in  peak  nodal  stress  predictions,  from  adjacent 
elements  at  the  hole  boundary,  is  less  than  0.5%. 
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4.2  Typical  stress  results  due  to  overload  followed  by  unload 

It  is  useful  to  consider  stress  distribution  results  due  to  overload  and  then  followed  by 
unloading  for  a  specific  case.  The  variation  of  plastic  zone  size  with  overload  level  is  then  also 
given. 

4.2.1  Typical  stress  result  due  to  overload  followed  by  unload 

The  distributions  of  the  stress  quantities  aVM  ,  (JYY ,  and  a  xx  resulting  from  a  compressive 
load  of  654  MPa  (equal  to  0.5  cry )  are  shown  in  Figure  7,  along  the  prospective  crack  path 
x=X / R.  For  this  overload,  the  plastic  zone  size  is  Dpz  =  0.3 R ,  as  indicated  by  the  plastic 
zone  size  in  Figure  7,  where: 

<JVM  is  the  von  Mises  stress; 

aYY  is  the  transverse,  v-direction  stress,  along  the  prospective  crack  path; 
dm  is  the  radial,  x-direction  stress,  along  the  prospective  crack  path. 


Nondimensional  distance  away  from  edge  of  hole,  x 


Figure  7:  Normalised  stress  components  along  prospective  crack  path  for  plate  with  remote  compressive 
overload 

After  the  overload  is  removed,  the  resultant  residual  stresses  along  the  prospective  crack  path 
where  y  =  0  are  as  shown  in  Figure  8.  There  is  an  inflection  point  in  <jyy  and  crVM  at  x  =  0.3. 

Also,  in  this  case  cryy  is  a  maximum  atx  =  0,  and  is  of  a  value  equal  to  0.54<ry  and  for  large x 
it  asymptotes  to  zero.  In  comparison  to  cryy,  the  u xx  stresses  are  low  and  hence  the 
magnitude  of  aVM  is  almost  the  same  as  cryy . 
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Figure  8:  Residual  stress  components  along  prospective  crack  path  for  compressive  overload  followed  by 
unload 

Consider  now  a  tensile  overload  of  654  MPa  followed  by  an  unload.  The  normalised  residual 
stress  (JYY  /  (7 Y  is  shown  in  Figure  9.  As  expected,  the  sign  of  the  value  is  the  opposite  to  that 
given  in  Figure  8,  but  the  magnitude  is  the  same. 


Figure  9:  Residual  transverse  stress  along  prospective  crack  path  for  large  D6ac  steel  plate  for  tensile 
overload  followed  by  unload 
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4.2.2  Variation  of  plastic  zone  size  with  remote  overload  stress 

Figure  10  gives  the  variation  of  the  nondimensional  plastic  zone  size  p=Dpz  /  R  with  the 
magnitude  of  the  compressive  overload.  This  is  a  useful  tool  for  determining  the  required 
initial  overload  stress  oini  for  any  particular  value  of  p  required.  This  plot  was  used  to  obtain 
values  of  <Jini  for  p  =  0.1,  0.2  and  0.3,  as  needed  in  later  sections. 


Nondimensional  compressive  overload  level,  3 a/ay 


Figure  10:  Depth  of  plastic  zone  versus  level  of  compressive  overload  for  a  round  hole  in  a 
unidirectional  stress  field 


4.3  Comparison  of  typical  stress  results  including  service  load 

Shown  in  Figure  11  is  the  nondimensional  transverse  stress  crYY  /  <jy  for  the  case  of  a  654  MPa 
tension  overload  followed  by  an  unload  and  then  application  of  a  200  MPa  service  load.  Also 
shown  are  the  results  for  the  654  MPa  compressive  overload  case.  As  expected,  the  tension 
overload  results  in  a  lowering  of  the  residual  stresses  aYY  /  crY  near  the  hole,  whereas  the 
compression  overload  results  in  higher  residual  stresses.  Note  also  that  the  edge  of  the  initial 
yield  zone  is  located  at  x  =  0.3.  As  expected,  for  large  x  the  two  solutions  converge. 
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0.0  0.5  1.0  1.5  2.0 


Nondimensional  distance  away  from  edge  of  hole,  x 

Figure  11:  Transverse  stress  for  overload  with  application  of 200  MPa  remote  load 

5.  Weight  function  method  for  cracked  hole  in  loaded 

plate 

This  section  initially  gives  the  definition  of  the  weight  function  approach  and  an  effective 
numerical  integration  technique  for  its  calculation.  Then  stresses  and  stress  intensity  factors 
for  the  cold-expanded-hole  case  are  given.  A  function  for  stress  at  the  edge  of  a  hole  due  to 
remote  loading  is  also  given  to  be  used  in  determining  the  effect  of  subsequent  remote  service 
load. 

5.1  Weight  function  definition 

Consider  an  initial  uncracked  hole  of  radius  R  in  a  loaded  plate  as  shown  in  Figure  12.  The 
Mode  I  stress  intensity  factor,  K  t ,  can  be  calculated  from  the  integral  of  the  product  of  the 
weight  function,  m(x,a) ,  and  the  transverse  distribution  of  stress  perpendicular  to  the 
prospective  crack  plane,  crvv(x) .  This  is  given  by  Bueckner  [1]  as: 

a 

K1  =  -JW  J  <jvv  (x)m(x,  a)dx  (2) 

o 
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where  W  is  a  characteristic  dimension,  and  for  cracks  emanating  from  a  hole  W  =  R . 


For  the  case  of  a  symmetric  crack  emanating  from  a  hole  in  a  large  plate,  it  has  been  shown  by 
Wu  and  Carlsson  [3]  that  this  weight  function  is  given  by: 


l  3  r 

m(x,a)  =  —j=  =2>,(o)  i- 
yJ27ra  ;=i  v 


a  y 


1-3/2 


(3) 


where  a  nondimensional  form  of  K  t  is  given  by: 

/=  Kr 


yjnaR 


(4) 


where: 

A  is  the  crack  length; 

a  is  the  crack  length  ratio,  A/R: 

x  is  the  normalised  distance  from  the  hole  edge,  X/R ; 

<7  iT  (x)  is  the  stress  in  the  uncracked  body  transverse  to  the  crack  path  at  a  distance  x  from 
the  hole  edge; 

J3j(a)  are  3  handbook  coefficient  values  (/  =  1  to  3)  taken  from  Wu  and  Carlsson  [3], 
Table  13.1,  with  N  =  2.  Here [ix (a)  =  2.0,  and  fiAa)  and  /!(«) are  functions  of  a, 
accurately  represented  by  ninth-order  polynomials. 


Figure  12:  Notation  and  geometry  for  analysis  of  overloaded  plate  containing  a  hole  and  local  detail  for 
calculation  of  stress  intensity  factor 

In  this  report  the  method  of  superposition  will  be  used  in  conjunction  with  the  weight 
function  method  for  determining  KI  due  to  combined  residual  stresses  and  stresses  from 
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subsequent  service  loading.  However,  in  some  cases,  negative  K j  will  be  determined.  The 
method  of  superposition  is  however  only  valid  when  the  crack  is  open.  The  limits  of  validity 
(of  crack  sizes  considered)  will  be  established  on  a  case-by-case  basis.  Here,  the  combined 
stresses  along  the  crack  path  determine  if  the  crack  is  fully  open  or  if  the  crack  is  partially 
closed.  For  this  formulation  of  the  weight  function  method  it  is  known  that  partial  closure  of 
the  crack  will  introduce  some  error  in  the  calculation  of  KI  depending  also  on  the  crack  size. 


5.2  Numerical  integration  method 

In  order  to  calculate  the  stress  intensity  factor  as  a  function  of  crack  length,  the  numerical 
method  by  Moftakhar  and  Glinka  [27]  is  used.  Reference  [27]  has  proposed  that  the  integral 
from  Equation  (2)  can  be  calculated  from  Equation  (5)  below: 

K,=Jv£JSiary(X)  (5) 

;=1 

This  represents  the  product  of  the  area,  Si  ,  under  the  curve  m(x,a )  and  the  value  of  the 
function  cr  vr  ( X ) .  And  this  corresponds  to  the  coordinate  X  of  the  centroid  of  the  area  St  as 
shown  in  Figure  12.  In  [27]  this  procedure  is  reversed  as: 

K,  =  yjw  ^S*m(X*  ,a)  (6) 

/=i 

This  procedure  is  implemented  by  calculating  the  areas  S/  under  the  stress  function 
avv  (x)  and  the  coordinates  X*  of  their  centroids.  This  leads  to: 

S*  =![>„(*,)  +  °yy (*m)]  (*<  “ *m)  (7) 

where  xt  is  the  x  coordinate  of  increment  i,  and  xi_l  the  x  coordinate  of  the  previous 
increment  and  X*  is  given  by: 


3  [^(x,)+CT,(xh)]  () 

The  weight  function  m(X* ,  a)  is  given  by  Equation  (3),  and  the  numerical  procedure  to  solve 
for  K f  has  been  implemented  in  a  Fortran  program  (see  Appendix  A).  Also  the  procedure 
adopted  here  uses  strips  of  width  dx  =  a  / 1 000  and  the  integration  is  taken  up  to  a  .  Since 
m(x,  a)  is  never  evaluated  any  closer  than  dx  1 2  from  x  =  a  the  singularity  in  m(x,  a)  at 
x  =  a  in  Equation  (3)  is  avoided. 
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It  should  be  noted  that,  in  order  to  represent  the  transverse  stresses  in  the  computer  program, 
a  4th-order  polynomial  was  used  in  the  initial  yield  zone  and  an  exponential  decay  was  used 
outside  this  region  (see  Figure  11). 


5.3  Comparison  for  remote  loading  case 

It  is  useful  to  confirm  the  accuracy  of  the  numerically  determined  value  of  K f ,  by  comparison 
with  a  known  analytical  solution,  for  the  remote-loading  case  that  is  shown  in  Figure  3.  For 
the  numerical  solution  the  equation  for  the  stress,  aYYre ,  away  from  the  edge  of  a  hole  in  an 
infinite  plate  for  elastic  conditions,  is  used.  This  is  given  by  Timoshenko  and  Goodier  [28]: 


cr 


YYre 


=  cr 


.  0.5 

1h - 7 

(1  +  x)2 


^  1.5 

+  (1  +  *)4 


(9) 


where  x  is  the  nondimensional  distance  ahead  of  the  hole  and  <jre  is  the  remote  stress. 

Hence,  the  normalised  stress  intensity  factor,/!  is  determined  numerically  using  Equations  (2), 
(3),  (4)  and  (9),  and  the  results  are  given  in  Figure  13.  The  corresponding  analytical  solution  by 
Wu  and  Carlsson  [3]  is  also  given,  showing  very  good  agreement. 


Figure  13:  Comparison  of  function  ffor  a  cracked  hole  in  loaded  plate  versus  crack  length  ratio 
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5.4  Stresses  and  stress  intensity  factor  for  cold-expanded  hole  (no  remote 
load) 

A  state  of  plane  strain  is  assumed  in  all  analyses  involving  the  cold-expanded  hole  unless 
noted  otherwise  in  order  to  compare  with  existing  analytical  solutions.  It  is  also  important  to 
note  that  an  elastic -perfectly  plastic  yield  condition  is  assumed  in  the  stress  response 
analysis1.  In  order  to  determine  the  nondimensional  initial  plastic  yield  zone  size, 
p  =  Dpz  /  R ,  corresponding  to  the  interference  level  A ,  the  following  equation  is  used  (see 
Equation  (64)  of  Jost  [29]): 


A=  <7y 


Ej 3 


(l  +  v)(l  +  /?)2 


f  pR  +  R  ^ 

2  \ 

l  +  (l-2v) 

V 

J 

+(1  +  vm)(1-2  vm)- 


^  M 


ln(l  +  p)  + 1  - 


pR  +  R 


W 


(10) 


where  vM  and  v  are  the  Poisson's  ratio  for  the  mandrel  and  plate  respectively,  and  EM  is  the 
modulus  of  elasticity  for  the  mandrel. 

5.4.1  Removal  of  mandrel  resulting  in  residual  stresses  for  case  of  no  re-yielding 

The  geometry  and  idealised  yielding  case  considered  here  is  shown  in  Figure  14.  It  is  assumed 
there  is  no  re-yielding  at  the  hole  edge  when  the  mandrel  is  removed. 


1  Hence,  some  small  discrepancy  is  expected  in  the  stress  prediction  as  compared  to  any  analysis  using 
the  strain  hardening  assumption  based  on  the  material  properties  in  Table  1. 
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Figure  14:  Geometry  and  idealised  yielding  cases  for  a  cold-worked  hole  in  a  remotely-loaded  finite- 
width  plate  with  no  re-yielding  following  removal  of  mandrel 


Rich  and  Impellizzeri  [30]  have  derived  the  residual  stress  distributions  around  a  hole  in 
which  a  mandrel  has  been  inserted  to  cause  plastic  deformation,  and  then  removed.  This 
induces  a  self -equilibrating  residual  stress  field.  For  convenience,  the  derivation  used  here  is 
given  by  Jost  [29] .  The  residual  tangential  stress,  aT  along  the  prospective  crack  path  v  =  0,  for 
an  elastic-perfectly  plastic  material  and  plane  strain  condition  for  x<  p  is  given  by: 


For  p<  x  <w: 


crvv(x)=crr 


[ji+pV  r/g+pflV  Jr  i  v 
|U+*J  I  w  J  L^1+Xy 


(12) 


From  Jost  [29]  and  converting  to  the  present  notation,  the  residual  radial  stress  for  x  <  p  is 
given  by: 
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cr . 


^xx(x)  =  ctr=-j^  pin 


^1  +  x^ 
\l  +  Pd 


-1  + 


R  +  pR 

V  w  J 


+  a 


Vl  +  xy 


R 

ywj 


1- 


Vwy 


(13) 


while  for  x  >  p  the  radial  stress  is  given  by: 


<v,M  =  ^=4N1+pT ''*+aR 


Vl  +  xy 


+  or 


y  w  j 


vl  +  xy 


/? 

W 


1- 


W 


(14) 


where  a  is  given  by: 


«  =  21n(l  +  p)  +  \- 


R(\  +  P) 

y  w  j 


(15) 


Another  Fortran  program  has  been  written  to  compute  f  the  nondimensional  stress  intensity 
factor,  using  weight  functions  from  Equations  (3)  and  (6)-(8)  and  Equations  (4)  and  (11)-(12). 
As  a  check  on  the  computer  code  accuracy,  results  for /obtained  using  the  present  method,  as 
R/w  approaches  zero,  are  shown  in  Figure  15,  together  with  the  corresponding  results  from 
Wu  and  Carlsson  [3].  Here,  good  agreement  is  seen. 


Figure  15:  Comparison  off  for  a  cold-worked  hole  versus  crack  length  ratio 
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5.4.2  Removal  of  mandrel  resulting  in  re-yielding 

Here,  we  consider  the  case  where  during  removal  of  the  mandrel  re-yielding  occurs  adjacent 
to  the  hole,  as  idealised  in  Figure  16.  A  peak  transverse  stress  exists  at  x  =  Rry  -  R  . 

It  can  be  shown  from  Jost  [29]  that  the  re-yield  radius  for  an  infinite  plate  is  given  by: 

RRY=^R\\  +  p)  e-°'25  (16) 

The  re-yield  radius  for  a  finite-width  plate  is  given  by: 


f 

-0.25  1h 

Rry=']r1({  +  P)  e 

rrt  )  f  w  J  y 

(17) 


This  allows  for  an  iterative  solution  for  Rry  .  The  nondimensional  re-yield  radius,  77 ,  is  given 
by: 


r/  =  (RRY-R)/R  (18) 

The  following  expressions  for  the  tangential  stress  are  all  derived  from  Jost  [29]  and  apply 
along  the  prospective  crack  path  where  y  =  0: 


a  =  A/R 
p  =  Dpz/R 
>1  =  (Rry~  R)/R 

4 -  2w  - ► 


Figure  16:  Geometry  and  idealised  yielding  cases  for  a  cold-worked  hole  in  a  remotely-loaded  finite- 
width  plate  with  re-yielding  following  removal  of  mandrel 
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For  0  <  x  <  t]  we  have: 


vjx)  =aT  =  - 


2a 


-p-(l  +  ln(l  +  x)) 
v  3 


(19) 


For  rj  <  x  <  p  we  have: 


cr 


cryy(x)  =  (Jr  =  ^p21n 


^1  +  x^ 
\l  +  Pj 


+  1  + 


RQ  +  P)' 

V  w  ) 


-2< 


V 


Vl  +  Xy 


+ 


Rrj 

) 


(20) 


And  for  p<x<w  we  have: 


..  cr  iM+pV  fi?(l  +  p) 

<7w(x)  =  <77.=-^-|  - -  +' 


VI 


-2< 


y  w  j 


V 


V  1  +  X  y 


+ 


Rr/ 

ywj 


(21) 


These  equations  are  sufficient  to  generate  the  entire  stress  solution  for  0  <  x  <  w  -  R  . 

5.4.3  Reaming  after  re-yielding 

Although  they  are  not  used  in  the  present  work,  for  completeness  the  stress  equations 
corresponding  to  reaming  after  re-yielding  are  also  given  here.  Usually  after  cold  working,  the 
hole  will  be  reamed  out,  increasing  the  size  of  the  hole  (to  create  a  smooth  undamaged 
surface).  The  following  equations  from  Jost  [29]  are  used.  As  before,  p  is  determined  from 
Equations  (1)  and  (10).  Firstly,  due  to  the  cold  expansion  we  have  the  real  radius  of  the  hole: 

Rreal={  l  +  QR  (22) 


Due  to  the  reaming  radius,  R RKAM  ,  we  have  the  increment  for  the  change  of  hole  radius: 

~  Rream  ~  Rreal  (23) 

From  Equation  (35)  of  Jost  [29]  we  have  the  following  equation  to  determine  the  new  re-yield 
radius,  Rry  : 


In ((R  +  da)  /  R)  =  2 In (R*ry)  +  (Rry  /  w)2  - (R*ry  /  wf  (24) 

The  nondimensional  re-yield  radius  is  given  by: 

r/  =  (R*RY-(R  +  Sa))/(R  +  Sa)  (25) 

The  following  stresses  apply  along  the  prospective  crack  path  where  v  =  0: 

For  rj  <  x  <  p  we  have: 
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2cr 


<7yy  (*)  =<7T=-  -J=~  (1  +  ln(l  +  X)  ) 


(26) 


For  rj  <  x  <  p  we  have: 


cr . 


ayy(x)=aT  =-j=]2\n 


< 1  +  x A 

v1  +  Py 


+  1  + 


(R  +  Ja)(l  +  p) 


w 


—  2< 


V 


v  l  +  xy 


+ 


( R  +  Sa)rj 


w 


(27) 


And  for  p  <  x  <w  we  have: 


cr 


%<+)  =  aT=-^ 


1  +  pY  |  f (/?  +  ffa)(l  +  /?) 

V  1  +  x  J  ^  w 


-2< 


7 


vl  +  xy 


+ 


(i?  +  ^a)^ 

y  w  j 


(28) 


These  equations  are  sufficient  to  generate  the  entire  stress  solution  for  0  <  x  <  vv-  /? .  It  has 
been  found  that  in  practice  very  little  difference  exists  between  the  unreamed  re-yielded 
solution  and  the  reamed  re-yielded  solution. 


5.5  Elastic  stress  for  finite-width  plate 

The  equations  above  for  the  tangential  stress  contain  a  finite-width  dimension  w ,  as  defined 
in  Figure  1.  However,  to  combine  with  an  elastic  remote  stress  also  requires  a  finite- width 
capability.  The  following  equation,  for  the  elastic  stress  in  the  neighbourhood  of  a  circular 
hole  in  a  finite-width  plate,  taken  from  Howland  [31],  is  used: 


°YYre  W  =  (T ,  J  1  +  2«t0  + 


d  n 


2X; 


n(2n  + 1  )d2n 
\  2/1+2 


( R(l  +  x)/wY  tURQ  +  x)/w) 


+ 


(n  ~  1)(2»  -  l)e2„ 
(R(l  +  x)/w)2'1 


+  n(2n  - 1  )l2n  (i?(l  +  x)  /  w)1"  2  +  (n  + 1)(2 n  + 1  )m2n  (i?(l  +  x)  /  w)2n  >  (29) 


where  mn,  dn,  en ,  ln  are  coefficients  defined  in  Howland  [31],  <Jre  is  the  remote  stress,  and  the 
normalised  distance  from  the  hole  edge  is  given  by  x  =  X/R.  The  only  limitation  is: 
l&YYre  +  GYYrS  \  -  ° y '  where  (J YYn,  is  obtained  from  Equation  (29)  above  and  <JYYrs  is  the  residual 
stress  defined  in  section  5.4.1  or  5.4.2  or  5.4.3.  As  a  check,  cr^  /  <rre  is  plotted  against 
R(\  +  x)/  w  for  different  R/ w  ratios,  as  shown  in  Figure  17.  Data  points  taken  from  [29]  are 
also  shown  in  Figure  17,  and  these  show  good  agreement  with  the  curves.  This  validates  this 
section  of  the  computer  code  developed  to  implement  the  analysis. 
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Figure  1 7:  Transverse  stress  <7m.e  /  a re  versus  R(l  +  x)/w  for  differen t  R/w,  using  analysis  given  in 
[31] 


6.  Cold-expanded  hole  subject  to  remote  loading  - 
weight  function  method 

In  this  section,  three  cases  for  cold-expanded  holes  are  considered  as  follows:  (i)  hole  in  a 
generic  large  D6ac  steel  plate;  (ii)  hole  in  a  generic  large  aluminium  alloy  plate  including 
comparison  with  published  results;  and  (iii)  application  to  a  hole  in  C-130  wing  skin  location, 
with  a  finite-width  geometry,  with  associated  checks  on  closure  effects.  All  cases  assume 
conditions  of  plane  strain. 


6.1  Application  to  a  hole  in  a  large  D6ac  steel  plate 

For  a  range  of  initial  plastic  yield  zone  sizes  using  Equations  (6)-(8),  the  tangential  stress 
variation  away  from  a  cold-worked  hole  has  been  determined  and  is  shown  in  Figure  18.  For 
completeness,  the  radial  stresses  are  also  given  in  Figure  19  for  various  nondimensional 
plastic  zone  sizes,  p  .  In  order  to  determine  the  range  of  validity  of  the  solution  for  the  stress 
intensity  factor,  the  sum  of  the  residual  stresses  and  those  due  to  subsequent  remote  loading 
stress  are  given  in  Figure  20,  along  the  prospective  crack  path.  The  remote  stress  is  200  MPa, 
and  it  is  seen  that  for  these  results,  for  p  =  0.1  and  0.2,  that  the  crack  is  fully  open  (i.e. 
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( 7yy  /  <Jv  >  0).  However,  for  p  =  0.3  there  is  crack  closure  between  0  <  a  <  0.084,  and  K t ,  as 

calculated  on  the  basis  of  a  fully-open  crack,  is  not  valid.  With  a  >  0.084,  the  crack  is  partially 
open  and  there  is  some  error  in  ,  which  reduces  for  increasing  a. 

Using  the  tangential  residual  stress  results,  a  Fortran  program  has  been  written  using  weight 
functions  to  compute  the  residual  stress  intensity  factor  caused  by  cold  working  of  the  hole  by 
a  mandrel.  These  results  are  shown  in  Figure  21  and  indicate  a  maximum  value  of  -763 
MPa  mm0-5  at  a  distance  of  0.6  mm  for  p  =  0.3.  For  smaller  initial  plastic  zone  sizes,  the 

maximum  value  of  K t  is  reduced  and  the  peak  value  occurs  at  a  smaller  crack  length  ratio. 
As  shown  in  Equations  (10)-(12),  the  smaller  the  initial  plastic  zone  size  the  smaller  the 
interference.  Similar  curves  have  been  obtained  by  Moreira  et  al  [20]  for  plane  stress 
conditions  also  using  a  weight  function  approach,  but  for  an  aluminium  alloy  plate  instead  of 
D6ac  steel. 


Figure  18:  Transverse  residual  stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate  due  to 
cold  working  followed  by  removal  of  mandrel 
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0.0  0.5  1.0  1.5  2.0 

Nondimensional  distance  away  from  edge  of  hole,  x 


Figure  19:  Radial  residual  stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate  due  to  cold 
working  followed  by  removal  of  mandrel 


Nondimensional  distance  away  from  edge  of  hole,  x 


Figure  20:  Transverse  stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate  due  to  cold 
working  and  effect  of  subsequent  remote  stress  of 200  MPa 
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Nondimensional  crack  length  ratio,  a 


Figure  21:  Stress  intensity  factor  results  for  a  cold-worked  hole  in  a  large  D6ac  steel  plate  from  residual 
stress 


To  investigate  the  beneficial  effects  of  cold  working  under  realistic  conditions,  a  combined 
loading,  using  the  method  of  superposition,  has  been  considered,  with  a  remote  tensile  stress 
of  200  MPa  and  initial  plastic  yield  zone  sizes  of  p  =  0.1, 0.2  and  0.3.  As  for  previous  cases,  the 
geometry  is  for  a  5  mm  radius  hole.  These  results  are  shown  in  Figure  22,  together  with  the 
results  for  the  case  of  remote  stress  alone.  Clearly,  the  cold  working  is  beneficial,  both  within 
the  initial  yield  zone  and  extending  some  distance  outside  the  initial  yield  zone,  with  the  best 
results  being  obtained  with  the  largest  initial  plastic  yield  zone. 
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0.0  0.5  1.0  1.5  2.0 

Nondimensional  crack  length  ratio,  a 


Figure  22:  Stress  intensity  factor  results  for  a  cold-worked  hole  in  a  large  D6ac  steel  plate  with 
subsequent  remote  tensile  stress  of 200  MPa 


6.2  Comparison  with  previous  work  for  a  large  aluminium  alloy  plate 

For  the  cold-worked  hole,  the  most  relevant  published  data  that  have  been  found  are  those  of 
Moreira  et  al  [20].  However,  they  correspond  to  an  infinite  aluminium  plate  under  plane  stress 
conditions.  Since  the  present  work  uses  conditions  of  plane  strain,  an  approximate 
comparison  was  made  by  adjusting  the  interference  in  the  present  analysis  in  order  to  achieve 
the  same  residual  stress  as  shown  in  Figure  23.  Here,  it  was  found  that  the  use  of  X  =  1.15% 
for  plane  strain  closely  agrees  with  A,  =  2.00%  for  the  plane  stress  case  analysed  in  Moreira  et 
al  [20],  The  geometry  used  is  presented  in  Table  2,  and  the  equations  used  were  Equations  (19) 
and  (21)-(23).  The  definition  of  interference  used  in  Moreira  et  al  [20]  is  the  same  as  that  in 
Equation  (1).  Also,  from  a  private  communication  with  the  corresponding  author,  it  was 
established  that  the  von  Mises  yield  criterion  was  used.  Note  that,  in  Figure  23,  in  the  set  of 
results  shown  for  Moreira  et  al  [20],  the  value  for  x  =  0.  This  corresponds  to  a ^  =  crY=  285 
MPa,  which  confirms  that  conditions  of  plane  stress  are  being  used.  Also,  while  re-yielding 
occurs  in  Moreira  et  al  [20],  which  also  occurs  in  our  analysis,  there  is  disagreement  atx  =  0, 
since  in  plane  strain  an  =  2a y  /  V3  =  329  MPa.  The  corresponding  residual  stress  intensity 
factor  for  both  the  present  work  and  that  from  Moreira  et  al  [20]  is  shown  in  Figure  24  and 
involves  a  disagreement  of  7.5%.  The  reason  for  the  disagreement  is  not  clear.  However,  in  the 
weight  function  method  used,  the  radial  stresses  will  not  contribute  to  the  stress  intensity 
factor,  when  using  a  J-integral  method  within  FEA  a  contribution  is  possible.  The  typical 
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normalised  tangential  and  radial  stresses  for  a  cold-worked  hole  in  an  aluminium  alloy  plate 
are  shown  in  Figure  25.  This  corresponds  to  1%  interference.  From  this  it  can  be  seen  that 
compressive  radial  stresses  exist  along  the  prospective  crack  path. 

Table  2:  Parameters  for  calculation  of  interference  and  assumptions  made 


Quantity 

Present  work 

Moreira  et  al  [20] 

em 

207000  MPa 

Not  applicable1 

EP 

71400  MPa 

71400  MPa 

285  MPa 

285  MPa 

R 

10  mm 

10  mm 

b 

finite  plate,  w/  P=10 

finite  plate,  w/  P=10 

0.3 

Not  applicable1 

Vp 

0.3 

0.3 

p 

0.875 

0.875 

A 

0.0115  (1.15%) 

0.02  (2%) 

Plane  strain 

Yes 

No 

Plane  stress 

No 

Yes 

Perfectly  elastic -plastic 

Yes 

Yes 

Re-yield  at  hole 

Yes 

Yes 

Yield  criteria 

von  Mises 

von  Mises2 

1  Interference  was  applied  as  an  imposed  displacement. 

2  From  private  communication  with  corresponding  author  (the  yield  criterion  that 
was  used  was  not  specified  in  journal  paper). 
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Figure  23:  Residual  stresses  along  prospective  crack  path  due  to  cold  expansion  in  a  large  aluminium 
plate,  comparison  with  work  in  [20] 


Figure  24:  Stress  intensity  factor  due  to  residual  stress  in  a  cold-worked  large  aluminium  plate  for 
1.15%  interference;  comparison  with  work  in  [20] 
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Figure  25:  Plot  of  nondimensional  tangential  and  radial  stresses  away  from  a  cold-worked  hole  in  an 
aluminium  plate 

6.3  Application  of  cold-expanded  hole  to  C-130  wing 

DSTO  is  embarking  on  a  substantive  coupon  test  program  in  support  of  likely  life  extension 
modifications  of  the  C-130  transport  aircraft  at  fatigue  hotspots.  A  typical  location  where  cold 
expansion  is  being  considered  is  at  the  lower  wing  skin  holes,  as  shown  in  Figure  26  and 
Figure  27.  Therefore,  to  assist  in  the  interpretation  of  the  pending  fatigue  test  results,  here 
predicted  K ,  results  are  determined  for  the  coupon  geometry  being  tested.  This  coupon  is 
made  from  an  aluminium  alloy  whose  properties  are  given  in  Table  3.  A  coupon  of  total  width 
2 w  =  32  mm  and  with  a  central  hole  of  radius  R  =  3.0  mm  is  to  be  tested  ( w / R  =  5.33),  and  it  is 
assumed  to  be  representative  of  this  structure. 


Table  3:  Properties  of  C-130  lower  wing  skin  at  CW1 


Material 

Young's 

Modulus 

Poisson's 

Ratio 

Compression 
Yield  Stress 

Hole 

Radius 

7075-T7351 

71000  MPa 

0.33 

372  MPa 

3.0  mm 

Based  on  the  cold-expansion  equations  in  Section  5.4.2,  shown  in  Figure  28  is  a  plot  of 
nondimensional  tangential  stress,  <JVV  (x)  /  crv ,  versus  the  nondimensional  distance  x  away 

from  the  hole  for  three  interference  levels.  The  original  yielded  zone  is  also  evident  along  with 
the  re-yielded  zone.  The  general  trend  is  for  higher  peak  values  (both  positive  and  negative) 
to  occur  for  increasing  levels  of  interference.  Furthermore,  the  peak  values  occur  at  increasing 
distances  away  from  the  hole  for  increasing  interference.  For  this  finite- width  plate,  the 
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positive  values  of  <7  (pc)  I  <Jy  peakatx  =  2.3  corresponding  to  5%  interference,  where  the  peak 
positive  nondimensional  stress  is  <Jyy  /  <jv  =  0.53. 


Figure  26:  Location  of  cold-expanded  hole  on  lower  wing  skin  ofC-130 
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CW-I 


Figure  27:  Local  detail  ofC-130  hole  subject  to  cold  expansion,  looking  in  spanwise  direction 


Nondimensional  distance  away  from  edge  of  hole,  x 

Figure  28:  Residual  stresses  along  prospective  crack  path  due  to  cold  expansion  in  a  finite-width 
aluminium  alloy  plate  representing  C-130 fatigue  test  coupons  (w/R  =  5.33) 

For  comparison  purposes,  corresponding  analyses  for  a  large  plate  with  w/R  =  3333  have  also 
been  undertaken,  in  order  to  simulate  the  behaviour  of  an  infinite  plate.  The  tangential  stress 
results  are  given  in  Figure  29.  It  can  be  seen  that  the  positive  values  of  nondimensional  stress 
< Tn,  (x)  /  crv  peak  atx  =  2.55  for  5%  interference  level,  where  the  peak  positive  nondimensional 

stress  is  cr  v  /  <7r  =  0.38.  Lienee,  it  is  apparent  that  the  peak  positive  stresses  are  significantly 
lower  for  the  simulated  infinite  plate  compared  to  those  for  the  finite-width  plate.  On  the 
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other  hand,  the  magnitudes  of  the  peak  negative  stresses  for  the  finite- width  plate  were  found 
to  be  marginally  lower  in  comparison  to  those  for  the  large  plate.  The  extent  of  the  zone  of 
reversed  yielding  is  also  smaller  for  the  finite- width  plate  than  for  the  simulated  infinite  plate. 
Hence,  it  is  important  that  the  finite  width  be  considered  when  interpreting  fatigue  test 
results. 

Although  the  analysis  presented  here  is  based  upon  2D  plane  strain  conditions,  it  is  useful  to 
compare  residual  stresses  with  a  full  3D  FEA  performed  by  Houghton  and  Campbell  [32].  For 
an  aluminium  alloy  plate  with  the  same  interference  level  of  3%,  the  same  residual  stress 
trends  are  predicted.  As  expected,  at  the  mid-plane  of  the  3D  analysis  good  agreement 
occurred  with  the  (plane  strain)  2D  peak  values,  but  the  x  coordinate  of  the  3D  peak  result 
was  lower.  At  the  surface  of  the  3D  analysis,  the  peak  values  were  lower;  however,  the  x 
coordinates  of  the  peaks  were  in  agreement. 

Using  the  weight  function  technique,  the  residual  stress  intensity  factor  obtained  for  the  plate 
with  w/R  =  5.33  representing  the  C-130  test  coupons  has  been  plotted  in  Figure  30.  The 
greatest  negative  stress  intensity  factor  is  achieved  with  the  highest  level  of  interference  of  5%, 
while  the  lowest  is  achieved  with  an  interference  level  of  3%. 


Nondimensional  distance  away  from  edge  of  hole,  x 


Figure  29:  Residual  stresses  along  prospective  crack  path  due  to  cold  expansion  in  a  simulated  infinite 
aluminium  alloy  plate  with  w/R  =  3333 
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Figure  30:  Residual  stress  intensity  factor  due  to  cold  working  of  aluminium  alloy  C-130  test  coupons 
with  w/R  =  5.33 

Now  consider  a  remote  applied  stress  that  is  equal  to  the  maximum  spectrum  tensile  stress  of 
153  MPa.  The  resulting  stress  intensity  factor  due  to  this  load  is  superimposed  on  the  residual 
stress  intensity  factor  and  is  shown  in  Figure  31.  The  beneficial  effect  of  the  higher  levels  of 
interference  is  obvious.  The  cold  expansion  is  very  effective  in  keeping  the  crack  closed  over  a 
significant  crack  length  ratio  or  minimising  K;  as  compared  to  the  no  cold-expansion  case. 
For  this  remote  load,  based  on  the  stresses  along  the  prospective  crack  path.  Figure  32 
indicates  that  the  crack  length  is  closed  for  0  <  a  <  0.85  for  the  5%  interference  case.  For  the  4% 
and  3%  interference  cases,  the  corresponding  sizes  are  0  <  a  <  0.73  and  0  <  a  <  0.57, 
respectively.  The  magnitude  of  the  error  for  large  cracks  has  not  been  fully  investigated  at  this 
stage. 

By  comparing  Figure  31  to  Figure  22,  it  can  be  seen  that,  for  the  realistic  cases  considered,  the 
cold  expansion  in  aluminium  alloys  is  much  more  effective  than  in  D6ac  steel  in  minimising 
Kj  values,  due  to  the  much  larger  initial  plastic  yield  zone  sizes. 
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Figure  31:  Stress  intensity  factor  due  to  cold  working  of  aluminium  alloy  C-130  fatigue  test  coupons 
with  w/R  =  5.33,  and  subsequent  remote  load  of  153  MPa. 
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Figure  32:  Transverse  stresses  along  prospective  crack  path  for  finite-width  aluminium  alloy  plate  with 
w/R  =  5.33  due  to  cold  working  and  subsequent  remote  stress  of  153  MPa 

6.3.1  Investigation  of  weight  function  integrand 

In  order  to  investigate  crack  closure  effects,  it  is  useful  to  explore  in  more  detail  the  weight 
function  integrand  <jvv  (x)m(x,  a )  for  the  case  of  cold  working  with  a  subsequent  remote  stress 

of  153  MPa.  This  is  undertaken  for  the  finite-width  plate  with  w/R  =  5.33  and  three  different 
crack  lengths,  i.e.  0.855  mm,  1.140  mm  and  1.710  mm,  for  the  3%  interference  case.  From 
Figure  33,  it  is  seen  that,  for  increasing  crack  lengths  beyond  a  =  0.57,  the  positive  values  of 
a  (. x)m(x,a)  increase  quickly  and  significantly.  For  this  particular  case,  this  behaviour 

suggests  that  the  negative  contribution  from  avv  (x)m(x,  a)  in  the  range  to  0.03  <  x  <  0.57,  due 
to  the  negative  values  of  avv  (x)  in  that  region,  can  be  ignored  for  large  cracks. 

Hence,  in  order  to  better  account  for  partial  crack  closure  when  using  the  weight  function 
method,  with  the  said  partial  crack  closure  being  caused  by  large  negative  stresses,  and  as  a 
conservative  approximation,  the  value  of  the  combined  stresses  <r  (x)  is  here  set  to  zero 

whenever  cr  T(x)  <  0 .  The  results  for  K  f  obtained  on  this  basis  are  shown  in  Figure  34.  In 
comparison  to  Figure  31,  apart  from  the  elimination  of  negative  values  of  K; ,  the  positive 
values  of  KI  have  been  increased.  A  possible  better  method  would  be  to  implement  the 
formulation  of  [25]  to  the  present  problem. 
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Figure  33:  Evaluation  of  the  product  of  the  transverse  stress  distribution  and  the  weight  function, 
<jvv  (x)m(x,  a ) ,  as  a  function  of  distance  away  from  the  hole  for  various  crack  lengths 


6.3.2  Accuracy  check  for  finite-width  correction  approach 

All  prior  sections  (i.e.  4,  5,  and  6.1)  have  considered  remote  loading  applied  to  a  large  plate. 
The  high  accuracy  of  the  implemented  weight  function  approach,  at  least  for  elastic 
conditions,  was  demonstrated  in  Section  5.3.  Here,  it  is  useful  to  also  check  the  accuracy  of  the 
weight  function  method  when  applied  to  a  finite- width  plate,  for  the  remote-load-only  case  of 
153  MPa  tension,  under  elastic  conditions.  A  standard  solution  for  K1  can  be  found  in  Rooke 
and  Cartwright  [33,  Figure  105],  corresponding  to  two  cracks  at  a  circular  hole  in  a  rectangular 
sheet  (sheet  height  equal  to  twice  the  sheet  width)  subjected  to  a  uniform  uniaxial  tensile 
stress.  The  data  are  quoted  to  be  accurate  to  within  1%  and  are  comprised  of  three  curves, 
corresponding  to  values  of  w/R  of  2,  4,  and  oo,  which  were  created  from  tabulated  values 
published  by  Newman  [34]. 

Approximate  solutions  for  the  present  case  with  w/R  =  5.33  for  nondimensional  crack  length 
ratios  of  a  =  0.5, 1.0,  and  1.5  were  obtained  from  [33]  by  interpolation.  The  comparison  of  the 
results  thus  obtained  with  the  weight  function  computations  is  given  in  Figure  34,  where  it 
can  be  seen  that  there  is  good  agreement  for  the  small  crack  size  of  a  =  0.5.  However,  for  crack 
lengths  that  are  longer  than  a  =  0.5,  an  increasingly  large  discrepancy  is  evident,  indicating 
that  K  j  is  here  being  underestimated.  This  can  be  expected,  because  the  weight  function  used 
here  is  applicable  to  cracks  in  large  plates,  but  has  been  utilised  here  in  conjunction  with  the 
stresses  associated  with  a  finite-width  plate.  Hence,  the  error  due  to  this  approximation 
increases  as  the  crack  size  increases.  It  should  be  noted  that  this  error  would  apply  linearly  to 
each  of  the  interference  level  cases,  and  the  relative  trends  would  be  unchanged. 
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Figure  34:  Stress  intensity  factor  due  to  cold  working  of  aluminium  alloy  C-130  test  coupons  with 
w/R  =  5.33  and  subsequent  remote  load  of  153  MPa  and  where  a  (x)  =  0  for  arv  <  0 


7.  Tensile  overload  case  -  weight  function  approach 

Consider  the  case  of  a  tension  overload  on  a  large  D6ac  steel  plate  containing  a  5  mm  radius 
hole  for  the  case  of  plane  stress.  Overloads  of  0.5a Y,  0.45<ry  and  0.4cry  followed  by 
unloading  are  investigated.  These  load  levels  are  sufficient  to  result  in  initial  plastic  yield  zone 
sizes  of  p  =  0.3,  0.2  and  0.1  respectively,  as  indicated  in  Figure  10. 

7.1  Residual  stress  field  for  uncracked  case 

The  residual  stress  field  has  been  determined  for  a  tension  overload  followed  by  an  unload  for 
three  different  initial  plastic  yield  zone  sizes  using  the  FEA  method  from  Section  4.  The 
normalised  transverse  residual  stresses  along  the  prospective  crack  path  are  shown  in  Figure 
35,  and  indicate  that  a  residual  compressive  stress  is  left  within  the  initial  plastic  yield  zone. 
Note  that  no  re-yielding  has  occurred  on  the  hole  surface  during  unloading.  The 
corresponding  normalised  radial  stresses  are  plotted  in  Figure  36,  and  these  are  of 
significantly  lower  magnitude  than  the  transverse  stresses. 

For  the  prospective  crack  path,  the  total  stress  due  to  residual  stress  and  a  remote  stress  of  200 
MPa  is  shown  in  Figure  37.  For  the  cases  with  p  =  0.1  and  0.2,  the  crack  is  always  open  for  the 

combined  condition  ( avv  /  ay  >  0 ).  Also,  it  can  be  seen  in  Figure  37  that  for  p  =  0.3  there  will 

be  crack  closure  for  0  <  x  <  0.031,  because  aYV  /  av  <  0  in  this  region.  Flence,  Kir  which  is 
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calculated  assuming  a  fully-opened  crack,  is  not  valid  in  this  range.  For  cracks  of  length 
x  >  0.031  there  is  some  error,  which  reduces  with  increasing  crack  length.  In  order  to 
investigate  this  further,  FEA  has  been  used  in  Figure  38  to  plot  the  percentage  of  crack 
opening  versus  nondimensional  crack  length,  a,  for  the  case  of  a  remote  load  of  200  MPa  and 
p  =  0.3.  This  shows  that,  up  till  crack  lengths  of  a  =  0.07,  the  crack  is  fully  closed,  and  for  crack 
lengths  of  a  >  0.1  the  crack  is  fully  open.  This  means  that  the  weight  function  method  will  be 
accurate  for  predictions  of  stress  intensity  for  a  >  0.1,  and  for  a  <  0.07  the  value  of  K ,  will  be 
zero.  Furthermore,  it  has  been  calculated  using  FEA  that  for  remote  stresses  of  200  MPa  the 
crack  is  fully  open.  Of  less  interest  is  the  180  MPa  case  for  p  =  0.3,  which  is  shown  in  Figure 
38,  which  was  only  considered  in  order  to  investigate  the  sensitivity  of  crack  closure  to  remote 
stress  level. 


Figure  35:  Transverse  residual  stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate,  due  to 
remote  tensile  overload  followed  by  unload 
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Figure  36:  Residual  radial  stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate,  due  to 
remote  tensile  overload  followed  by  unload 


Figure  37:  Transverse  stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate;  total  stress  due 
to  subsequent  remote  load  of 200  MPa 
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Figure  38:  Length  of  crack  opening  versus  nondimensional  crack  length  ratio,  a,  for  remote  load  of  200 
MPa  and  180  MPa,  where  p  =  0.3 


Figure  39:  Stress  intensity  factor  for  a  large  D6ac  steel  plate  following  tension  overload  and  unload 
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7.2  Stress  intensity  factor  due  to  residual  stress  field 

Using  the  weight  function  technique,  from  Equation  (2),  the  residual  stress  intensity  factors 
are  determined  for  nondimensional  plastic  zone  sizes  corresponding  to  p  =  0.1,  0.2  and  0.3. 
The  residual  stress  intensity  versus  crack  length  ratio  is  given  in  Figure  39,  and  for  p  =  0.3  it 
shows  a  peak  negative  stress  intensity  of  569  MPa  mm0-5  at  a  crack  length  ratio  of  a  =  0.09. 
Smaller  initial  yield  zone  sizes  show  a  reduced  peak  negative  stress  intensity  factor  and  the 
peak  occurs  at  a  lower  crack  length  ratio. 

7.3  Stress  intensity  factor  due  to  overload  and  subsequent  remote  loading 

Stress  intensity  factor  results  are  given  here  for  the  combined  case  of  overload-induced 
residual  stress  with  the  addition  of  a  remote  tensile  load  of  200  MPa,  as  shown  in  Figure  40. 
Once  again,  nondimensional  plastic  zone  sizes  of  p  =  0.1,  0.2  and  0.3  have  been  considered. 
When  p  =  0.3,  the  edge  of  the  yielded  zone  extends  out  to  a  =  0.30.  For  comparison  purposes, 
the  results  for  the  case  of  the  remote  tensile  load  alone  are  also  given.  The  beneficial  effects  of 
the  residual  stresses  due  to  the  overloading  are  clearly  evident  when  the  crack  is  relatively 
small.  However,  once  the  crack  has  grown  through  the  initial  yield  zone  the  beneficial  effect  is 
eliminated.  As  expected,  greater  beneficial  effects  are  apparent  for  the  larger  initial  plastic 
zone  sizes.  As  a  check  on  the  accuracy  of  the  present  approach,  results  for  the  case  of  remote 
stress  only  are  compared  to  a  linear-elastic  analytical  solution  taken  from  Bowie  [35].  The 
comparison  is  given  in  Figure  40  where  very  good  agreement  is  shown  for  the  remote  load 
only  case. 


Nondimensional  crack  length  ratio,  a 


Figure  40:  Stress  intensity  factor  results  due  to  remote  overload  for  a  large  D6ac  steel  plate  with 
subsequent  remote  load  of 200  MPa 
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8.  Compression  overload  case  -  weight  function 

approach 

In  this  section,  we  consider  a  large  remotely-loaded  D6ac  steel  plate  with  a  central  hole  (see 
Figure  1).  Initially,  the  residual  stress  field  is  determined  by  FEA,  for  the  uncracked  case  for  2 
levels  of  initial  plastic  yield  zone  size,  due  to  a  remote  compressive  overload.  Then  the 
corresponding  stress  intensity  factor  results  are  given.  Finally,  the  stress  intensity  factor  is 
determined  with  a  subsequent  remote  tensile  service  load  of  200  MPa. 

8.1  Residual  stress  field  for  uncracked  case 

The  normalised  transverse  residual  stresses,  avv  /  cr, ,  along  the  prospective  crack  path  have 
been  computed  for  two  values  of  nondimensional  plastic  zone  size,  p  =  0.2  and  0.3,  and  are 
shown  in  Figure  41.  In  this  instance,  it  is  apparent  that  the  compression  overload  has  left 
residual  tensile  stresses  within  the  initial  plastic  yield  zone.  As  expected,  they  are  the  same 
magnitude  but  of  opposite  sign  to  those  given  in  Figure  35  for  the  tensile  overload  case. 

8.2  Stress  intensity  factor  due  to  residual  stress  field 

From  the  residual  stress  data  in  Figure  41,  Equation  (2)  has  been  used  to  generate  the  stress 
intensity  factor  as  a  function  of  nondimensional  crack  length  ratio,  a,  and  the  results  are 
shown  in  Figure  42.  Here  the  maximum  residual  value  of  is  559  MPa  mm° 5  for  p  =  0.3 
and  359  MPa  mm0-5  for  p  =  0.2. 


Figure  41:  Stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate,  residual  stresses  due  to 
compressive  overload  followed  by  unload 
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Figure  42:  Stress  intensity  factor  results  due  to  remote  compression  overload  for  a  large  D6ac  steel 
plate  followed  by  unload 


Figure  43:  Total  stresses  along  prospective  crack  path  for  a  large  D6ac  steel  plate  due  to  remote 
compression  overload  followed  by  unload  and  then  subsequent  remote  tensile  stress  of  200 
MPa 
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8.3  Stress  intensity  factor  due  to  overload  and  subsequent  remote  loading 

The  total  stress  due  to  the  overload-induced  residual  stress  and  a  remote  stress  of  200  MPa  is 
shown  in  Figure  43  for  the  uncracked  body.  As  expected,  it  can  be  seen  that  a  prospective 
crack  would  be  open  at  all  times  and  the  weight  function  method  will  be  valid  for  all  crack 
lengths.  The  corresponding  stress  intensity  factors  are  shown  in  Figure  44.  Flere  the  maximum 
Kj  of  1279  MPa  mm0-5  occurs  at  a  crack  length  ratio  of  a  =  0.14,  i.e.  0.7  mm,  for  an  initial  yield 
zone  size  of  p  =  0.3.  For  p  =  0.2,  the  maximum  Kl  is  1022  MPa  mm0-5  and  occurs  at  a  crack 
length  ratio  of  a  =  0.12,  i.e.  0.6  mm.  However,  when  the  growing  crack  extends  beyond  the 
initial  yield  zone,  K;  reduces  significantly  and  closely  follows  the  curve  for  the  case  of  no 
residual  tensile  stress  but  with  remote  loading  only. 


Non-dimensional  crack  length  ratio,  a 


Figure  44:  Stress  intensity  factor  results  due  to  remote  compression  overload  for  a  large  D6ac  steel 
plate  with  subsequent  remote  stress  of 200  MPa 


9.  Comparison  of  tension  overload  and  cold  expansion 

In  this  section,  initially  comparison  is  made  between  plane  stress  and  plane  strain  results 
obtained  for  the  tension  overload  case  for  a  large  D6ac  steel  plate.  Then  comparison  between 
cold  expansion  and  tension  overload  for  the  same  plastic  zone  size  is  made,  again  for  the  case 
of  a  large  D6ac  steel  plate.  Furthermore,  a  comparison  is  made  between  a  tension-overloaded 
and  cold-expanded  hole  in  a  large  aluminium  plate  for  the  case  of  plane  strain. 
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Figure  45:  Comparison  of  transverse  stresses  aYY  /  <jy  for  plane  stress  and  plane  straw  for  tension 
overload  followed  by  unload  in  a  large  D6ac  steel  plate 

9.1  Comparison  of  stress  intensity  factors  for  plane  stress  and  plane  strain 

Most  of  the  overload  cases  reported  up  to  this  stage  have  been  undertaken  assuming  plane 
stress  while  the  equations  for  cold  expansion  correspond  to  the  case  of  plane  strain.  Hence,  it 
is  useful  to  first  make  a  comparison  between  plane  stress  and  plane  strain  on  the  basis  of  the 
same  plastic  zone  size,  p ,  for  D6ac  steel,  both  with  a  tension  overload,  as  shown  in  Figure  45. 
It  can  be  seen  that  the  plane-strain  case  results  in  slightly  higher  positive  peak  values  and  also 
higher  negative  peak  values  at  the  same  p  value. 

The  computed  residual  stress  intensity  factor  due  to  the  overload  is  shown  in  Figure  46  under 
conditions  of  plane  stress  and  plane  strain  with  the  same  plastic  zone  size.  The  residual  stress 
intensity  factor  displays  a  considerably  greater  peak  negative  value  for  plane  strain,  being  - 
808  MPa  mm0-5,  compared  with  -571  MPa  mm0-5  for  the  case  of  plane  stress. 

Consider  now  the  stress  intensity  factor  subject  to  a  service  load  of  200  MPa,  as  shown  in 
Figure  47.  Except  for  very  small  crack  lengths,  the  disagreement  is  small.  However,  the 
overload  that  is  required  to  produce  the  specified  plastic  zone  size  is  significantly  different, 
being  a  =  0.5crv  for  plane  stress  and  cr  =  0.622cr .  for  plane  strain.  This  is  expected  because 
the  von  Mises  stress  yield  condition  depends  on  the  triaxial  stress  field. 
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Figure  46:  Comparison  of  KI  for  plane  stress  and  plane  strain  for  tension  overload  in  a  large  D6ac 
steel  plate 
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Figure  47:  For  tension  overload  a  comparison  of  KI  for  plane  stress  and  plane  strain  for  the  same 
initial  plastic  yield  zone  size  in  a  large  D6ac  steel  plate,  subject  to  service  load  of 200  MPa 
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Figure  48:  Comparison  of  stresses  along  prospective  crack  path  for  tension  overload  and  cold-worked 
hole  subsequent  to  remote  load  in  a  large  D6ac  steel  plate,  assuming  plane  strain 

9.2  Comparison  of  results  for  cold  expansion  and  tension  overload 

It  is  useful  to  compare  the  results  obtained  from  cold  expansion  with  those  from  tension 
overload,  since  they  both  result  in  residual  compressive  stresses  at  the  hole  edge.  However,  all 
the  overload  cases  presented  in  the  previous  section  are  for  plane  stress  conditions,  while  the 
equations  for  cold  expansion  were  derived  for  the  case  of  plane  strain. 

9.2.1  Steel  plate  example 

To  allow  comparison,  the  previous  FEA  is  repeated  for  the  tension-overload  case  but  this  time 
using  plane  strain  conditions  in  a  large  plate  with  a  hole  of  radius  R  =  5  mm.  The  amplitude  of 
the  remote  tension  load  is  chosen  so  that  the  plastic  zone  size  is  p  =  0.3.  The  comparison  of 
the  residual  transverse  stress  results  is  given  in  Figure  48  for  D6ac  steel.  This  shows  that  the 
tension-overload  case  has  a  much  higher  positive  peak  stress  than  the  cold-expansion  case. 
The  peak  negative  stress  value  is  also  larger  for  the  tension-overload  case.  The  residual  stress 
intensity  factor  is  shown  in  Figure  49  and  indicates  that  for  short  cracks  (i.e.  for  values  of  a  in 
the  range  0.05-0.13)  the  tension  overload  is  more  effective.  However,  the  cold-expanded  hole 
is  significantly  more  effective  away  from  the  edge  of  the  hole. 

These  differences  are  expected  to  be  due  to  the  fact  that  the  cold  working  affects  the  entire 
circumference  of  the  hole,  while  the  effect  of  the  tension  overload  is  more  localised.  Hence,  the 
tension  overload  is  only  expected  to  be  useful  for  relatively  small  cracks,  as  compared  to  the 
hole  size.  Some  further  minor  contribution  to  these  differences  can  also  be  expected  since,  for 
the  tension-overload  case,  the  material  response  has  some  post-yield  strain  hardening  and  the 
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plate  width  is  R/W=  10,  while  for  the  cold-expansion  case  R/W approaches  infinity,  and  the 
material  response  has  no  post-yield  hardening. 


Nondimensional  crack  length  ratio,  a 


Figure  49:  Comparison  of  residual  stress  intensity  factor  for  tension-overloaded  hole  and  cold-worked 
hole  in  a  large  D6ac  steel  plate 


Nondimensional  crack  length  ratio,  a 

Figure  50:  Comparison  of  stress  intensity  factors  for  a  large  D6ac  steel  plate  with  remote  tension 
overload  or  cold-worked  hole  and  the  same  subsequent  remote  load  of 200  MPa 
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A  comparison  of  results  between  a  tension-overloaded  hole  case  and  a  cold-expanded  hole 
case,  both  subjected  to  a  remote  tension  stress  of  200  MPa,  has  been  made  in  Figure  50  for 
D6ac  steel.  It  can  be  seen  that  cold  working  is  generally  more  effective  in  minimising  K {  once 
a  >  0.14,  and  has  an  influence  beyond  the  initial  plastic  yield  zone.  This  may  be  due  to  the  fact 
that  a  cold-worked  hole  has  plasticity  induced  all  the  way  around  the  hole,  while  the  tension 
overload  produces  localised  plasticity  only. 

9.2.2  Aluminium  alloy  plate  example 

To  further  explore  these  trends,  the  case  of  an  aluminium  alloy  plate  is  considered  for  plane 
strain  conditions.  In  Figure  51  a  comparison  is  made  between  a  tension-overloaded  (320  MPa) 
and  cold-expanded  hole  (2.02%)  in  aluminium  for  p  =  1.193.  In  both  cases  plane  strain  is 
assumed  and  a  subsequent  remote  tension  stress  of  153  MPa  is  then  applied.  The  cold- 
expanded  hole  is  significantly  more  effective  in  reducing  KI  than  the  overloaded  hole.  Note 
also  that  there  is  a  local  increase  in  the  overload  curve  at  a  crack  length  ratio  of  a  =  0.15;  this  is 
due  to  re-yielding  at  a  distance  corresponding  to  a  nondimensional  re-yield  radius  of  T) «  0.15. 


Figure  51:  Comparison  of  stress  intensity  factors  for  a  large  aluminium  alloy  plate  with  remote  tension 
overload  or  cold-worked  hole  and  the  same  subsequent  remote  load  of  153  MPa,  for  plane 
strain 
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10.  Conclusions 

A  weight  function  approach  has  been  used  to  determine  stress  intensity  factors  for  cracks  at 
holes  in  metallic  plates,  in  the  presence  of  either  compressive  or  tensile  residual  stress  fields, 
caused  by  three  types  of  overloads,  and  with  and  without  the  subsequent  application  of 
remote  loading.  The  overload  types  are  remote  tension,  remote  compression  and  hole  cold 
expansion.  Many  numerical  cases  are  considered  to  investigate  the  key  trends  in  stress 
intensity  factor. 

Some  points  to  note  are  as  follows: 

(i)  All  stress  results  in  the  corresponding  uncracked  body  are  presented  in 
nondimensional  terms  normalised  to  the  yield  stress,  so  there  is  some  transferability 
of  the  results. 

(ii)  Cold  expansion  of  both  steel  and  aluminium  plates  for  realistic  interference  levels 
was  investigated.  It  was  found  that,  for  the  same  plastic  zone  size,  cold  expansion  is 
more  effective  in  reducing  the  stress  intensity  factor  than  remote  tensile  overload. 

(iii)  For  realistic  cold-expansion  levels,  a  much  larger  plastic  zone  size  was  achieved  by 
cold  expansion  of  a  hole  in  an  aluminium  alloy  plate  as  compared  to  a  D6ac  steel 
plate,  resulting  in  an  increased  relative  benefit  in  reducing  the  stress  intensity 
factor. 

(iv)  For  all  overload  cases,  for  cracks  in  the  initial  yield  zone  a  greater  effect  on  stress 
intensity  factor  occurred  for  larger  magnitudes  of  the  overload. 

(v)  For  both  remote  overload  cases,  it  was  shown  that,  once  the  crack  length  was  the 
same  or  larger  than  the  initial  yield  zone  size,  the  stress  intensity  factor  was  the 
same  as  for  the  case  without  the  initial  overload. 

(vi)  Flowever,  for  the  cold-expanded  hole  case,  the  beneficial  reduction  in  stress 
intensity  factor  extended  (as  compared  to  the  tensile  overload)  for  some  distance 
outside  the  initial  yield  zone. 

(vii)  The  application  of  cold  expansion  to  an  aluminium  alloy  fatigue  test  coupon 
representing  a  lower  wing  skin  location  in  the  C-130  aircraft  was  considered.  For  a 
finite- width  plate,  the  compressive  residual  stresses  at  the  hole  edge  were  greater  in 
magnitude  than  for  a  large  plate. 

(viii)  Overall  the  present  results  have  been  benchmarked  favourably  against  other 
solutions  where  they  are  available. 

(ix)  For  cases  with  local  residual  compressive  stresses,  some  understanding  has  also 
been  gained  of  the  potential  inaccuracies  in  the  present  formulation  of  the  weight 
function  method  due  to  the  presence  of  regions  of  likely  crack  closure.  This  effect 
reduces  with  increasing  crack  length.  It  is  recommended  that  further  work  be 
undertaken  to  implement  the  method  of  Yu  and  Abel  [25]  to  the  present  problem. 
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Appendix  A:  List  of  computer  programs 

A  set  of  self-contained  Fortran  computer  programs  was  written  to  generate  the  data  in  the 
Figures,  and  these  programs  are  contained  on  a  Compact  Disc.  These  programs  were  run  on 
the  DSTO  server  Shyted  under  Fortran  77.  The  names  of  the  programs,  together  with  the 
Figure  numbers  in  whose  generation  they  participated,  are  presented  in  Table  Al.  The 
contents  of  the  files  on  Compact  Disc  are  described  in  Table  A2.  The  variables  used  in  the 
computer  programs  are  described  in  Table  A3. 


Table  Al:  Fortran  computer  programs  used  to  generate  data  in  the  Figures 


Program  Number 

Program  Name 

Name  Of  Main 
Output  File 

For  Generation  Of 
Data  In  Figure 

25 

wtfuncOO.f 

ansi 

13 

26 

wtfunc50.f 

ans 

15 

24 

wtfunll50.f 

what 

17 

36 

wtfungjlO.f 

new 

18 

37 

wtfungjll.f 

new 

19 

38 

wtfungjl2.f 

new 

20 

1 

wtfunll.f 

ansO 

21 

1 

wtfunll.f 

ansO 

22 

40 

wtfungjps3.f 

new 

23 

41 

wtfungjps4.f 

ans 

24 

39 

wtfungj.f 

new 

25 

33 

wtfungj05.f 

new 

28 

31 

wtfungjOl.f 

new 

29 

32 

wtfungj04.f 

ans 

30 

32 

wtfungj04.f 

ans 

31 

34 

wtfungj06.f 

new 

32 

22 

wtfunllO.f 

valid 

33 

23 

wtfunlll.f 

ans 

34 

16 

wtfunl03.f 

check 

37 

17 

wtfunl04.f 

check 

37 

18 

wtfunl05.f 

check 

37 

19 

wtfunl06.f 

ans 

39 

20 

wtfunl07.f 

ans 

39 

21 

wtfunl08.f 

ans 

39 

19 

wtfunl06.f 

ans 

40 

20 

wtfunl07.f 

ans 

40 

21 

wtfunl08.f 

ans 

40 

29 

wtfuncx.f 

ansi 

42 

2 

wtfun33.f 

ansO 

42 

14 

wtfunlOl.f 

checkl 

43 

15 

wtfunl02.f 

check2 

43 

30 

wtfuncy.f 

ans 

44 

27 

wtfunch.f 

ans 

44 

28 

wtfunci.f 

ans 

44 

6 

wtfun78.f 

check 

45 
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Table  A1  (continued):  Fortran  computer  programs  used  to  generate  data  in  the  Figures 


Program  Number 

Program  Name 

Name  Of  Main 
Output  File 

For  Generation  Of 
Data  In  Figure 

7 

wtfun79.f 

checkl 

45 

9 

wtfun81.f 

ans6 

46 

8 

wtfun80.f 

ans5 

46 

4 

wtfun74.f 

ans8 

47 

5 

wtfun75.f 

ansO 

47 

10 

wtfun82.f 

check 

48 

11 

wtfun83.f 

checkl 

48 

13 

wtfun85.f 

ansO 

49 

12 

wtfun84.f 

ansi 

49 

5 

wtfun75.f 

ansO 

50 

1 

wtfunll.f 

ansO 

50 

3 

wtfun73.f 

ans 

51 

42 

wtfungjps5.f 

ans 

51 

35 

wtfungj07.f 

ans 

— 
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Table  A2:  Description  of  files  on  Compact  Disc. 


Program 

Description 

Figures 

wtfuncOO.f 

Two  symmetric  cracks  emanating  from  a  hole  subject 
to  remote  tensile  stress.  (Validation  of  weight  function 
technique.) 

13 

wtfunc50.f 

Cold-worked  hole.  (Validation  of  Wu  and  Carlsson 
derivation.) 

15 

wtfunll50.f 

Howland  plot,  using  Equation  (29) 2,  obtained  from 
Howland  [31]. 

17 

wtfungjlO.f 

Cold-worked  hole  and  subsequent  remote  loading. 
Uses  Equation  (9),  obtained  from  Timoshenko  and 
Goodier  [28], 

18, 19,  20,  21,  22 

wtfungj04.f 

Cold-worked  hole  and  subsequent  remote  loading. 
For  finite-width  plate.  Uses  Equation  (29),  obtained 
from  Howland  [31].  Also  has  re-yielding  or  no  re¬ 
yielding  case. 

28,  29,  30,  31,32 

wtfunlll.f 

As  above,  except  that  =  0  for  cr  <  0  for  re-yielding 

analysis  only. 

34 

wtfungj07.f 

Ream  after  re-yield.  Cold-worked  hole  and 
subsequent  remote  loading.  For  finite-width  plate. 
Uses  Equation  (29),  obtained  from  Howland  [31]. 

wtfunl03.f 

Tension  overload  case.  Results  from  FEA  used  to 
generate  fourth-order  polynomial  and  exponential 
decay  curve  fit  for  p  =  0.3.  Uses  Equation  (9), 
obtained  from  Timoshenko  and  Goodier  [28]. 

37,  39,  40 

wtfunl04.f 

Tension  overload  case.  Results  from  FEA  used  to 
generate  fourth-order  polynomial  and  exponential 
decay  curve  fit  for  p  =  0.2.  Uses  Equation  (9), 
obtained  from  Timoshenko  and  Goodier  [28]. 

37,  39,  40 

wtfunl05.f 

Tension  overload  case.  Results  from  FEA  used  to 
generate  fourth-order  polynomial  and  exponential 
decay  curve  fit  for  p  =  0.1.  Uses  Equation  (9), 
obtained  from  Timoshenko  and  Goodier  [28]. 

37,  39,  40 

wtfuncx.f 

Compression  overload  case.  Results  from  FEA  used  to 
generate  third-order  polynomial  and  exponential 
decay  curve  fit  for  p  =  0.3.  Uses  Equation  (9), 
obtained  from  Timoshenko  and  Goodier  [28]. 

42,  43,  44 

wtfun33.f 

Compression  overload  case.  Results  from  FEA  used  to 
generate  third-order  polynomial  and  exponential 
decay  curve  fit  for  p  =  0.2.  Uses  Equation  (9), 
obtained  from  Timoshenko  and  Goodier  [28]. 

42,  43,  44 

2  Refers  to  equation  number  as  identified  in  present  report. 
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Table  A3:  Variables  used  in  computer  programs 


Variable 

Description 

Sig 

Yield  stress 

Sigl 

Remote  applied  stress 

E 

Young's  modulus  of  plate 

Em 

Young's  modulus  of  mandrel 

V 

Poisson's  ratio  for  plate 

vm 

Poisson's  ratio  for  mandrel 

R 

Radius  of  hole 

rho 

*Nondimensional  plastic  yield  zone  (p  in  report) 

RH 

*Length  of  yielded  zone  (Dp/  in  report) 

bb 

Finite  width  (w  in  report) 

Lamda 

Interference  ( A  in  report) 

*  Is  calculated  in  cold  expansion  programs,  but  must  be  set  in  overload  program. 
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For  life  assessment  of  airframe  components,  it  is  important  to  understand  how  the  residual  stresses  due  to  material  yielding  affect  crack 
initiation  or  crack  growth  at  holes.  A  key  step  in  understanding  such  effects  is  the  quantification  of  the  Mode  1  stress  intensity  factor  for 
cracks  that  may  initiate  subsequent  to  overloading.  Here,  a  two-dimensional  weight  function  approach  is  used  to  determine  stress  intensity 
factors  for  cracks  in  either  tensile  or  compressive  stress  fields  due  to  one  of  three  mechanisms:  remote  tension  overload,  remote  compression 
overload  or  hole  cold  expansion.  The  effect  of  subsequent  remote  loading  is  also  considered.  The  key  input  is  the  stress  distribution  in  the 
corresponding  uncracked  body  along  the  prospective  crack  path.  The  key  trends  are  investigated  through  many  numerical  examples  for 
symmetrically-cracked  holes  in  large  steel  and /  or  aluminium  alloy  plates.  Cold  expansion  of  finite-width  plates  representative  of  C-130  wing 
skin  locations  is  also  studied.  For  both  remote  overload  cases,  it  is  shown  that,  once  the  crack  length  is  the  same  or  larger  than  the  initial  yield 
zone,  the  stress  intensity  factors  are  the  same  as  for  the  case  without  the  initial  overload.  However,  for  cold-expanded  holes,  the  beneficial 
reduction  in  stress  intensity  factor  extends  beyond  the  initial  yield  zone.  Hence,  the  present  work  provides  a  greater  depth  of  understanding 
of  how  typical  residual  stresses  can  affect  key  inputs  into  airframe  life  assessment. 
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